Lifetime and decay of unstable particles in strong gravitational fields 
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We consider here the decay of unstable particles in geodesic circular motion around compact 
objects. For the neutron, in particular, strong and weak decay are calculated by means of a semi- 
classical approach. Noticeable effects are expected to occur as one approaches the photonic circular 
orbit of realistic black- holes. We argue that, in such a limit, the quasi-thermal spectrum inherent 
to extremely relativistic observers in circular motion plays a role similar to the Unruh radiation for 
uniformly accelerated observers. 



PACS numbers: 95.30.Cq, 14.20.Dh 



I. INTRODUCTION 

The possible decay of inertially stable particles due to 
strong gravitational fields has been considered recently. 
In particular, the proton decay, by weak and strong in- 
teractions, in uniformly accelerated trajectories [1 0, Q 
and in circular motion around compact objects[4|, has 
been considered in great detail. The astrophysical im- 
plications of these results are now under investigation. 
The possible decay of accelerated protons, however, is 
not a new issue. It can be traced back to the works of 
Ginzburg and Zharkov[5j in the sixties, where processes 
of the type p + — —> n°Tr + were considered. At the same 
time, ZharkovQ investigated the weak and strong de- 
cay of protons accelerated by an external electromagnetic 
field. (See for a review.) Clearly, none of these pro- 
cesses would be allowed in the absence of external forces. 
We notice, however, that there are subtle differences be- 
tween processes involving uncharged particles where the 
accelerations have gravitational and electromagnetic ori- 
gins, see [3] for further details. 

In this paper, we consider the decay and the lifetime 
of unstable particles in geodesic circular motion around 
spherically symmetrical compact objects. We evaluate, 
in particular, decay rates and lifetime for neutrons in 
relativistic circular motion according to the semiclassical 
approach introduced in [J]. Both the weak 
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channels are considered. Our results are compared to 
several ones obtained previously in the literature for the 
related processes p + — — > n a e + v e and p + — — ► n°n + . As 
we will see, for realistic black-holes, noticeable effects are 
expected to occur for circular geodesies close to the pho- 
tonic orbit r = iGM/rc 2 . Observers in these (unstable) 
circular orbits are necessarily in extremely relativistic 



motion (v 



w ), and it is well known that they indeed 



realize the inertial vacuum as a quasi-thermal distribu- 
tion ofoparticles characterized by a temperature T in the 
range ; 



ha , „ ha 
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where a stands for the effective Minkowskian centripetal 
acceleration for relativistic circular orbits. Our results 
suggest that the temperature ([3]) have in the present case 
the same central role played by Unruh temperature flfj| 
in the analysis of uniformly accelerated particles as seen 
from Rindler observers @- (See [§| for a recent review.) 
Similar conclusions hold also for other unstable particles. 



II. CIRCULAR GEODESICS AROUND 
COMPACT OBJECTS 

The line element corresponding to a spherically sym- 
metrical object of mass M is given by the Schwarzschild 
metric 

ds 2 = - (1 - 2M/r) dt 2 + (1 - 2M/r)~ X dr 2 + r 2 d£ 2 , 

(4) 

where dY? = d6 2 + sin 9 dtp 2 . Natural unities are adopted 
hereafter. In this coordinate system [HI]. a particle of 
mass m in a circular timelike geodesic at radius r on the 
equatorial plane have energy per mass ratio given by 



S/m = (1 - 2M/r)i = (1 - 2M/r)/ a/I - 3M/r, (5) 

with the dot standing for s-derivative. Its angular mo- 
mentum L = r 2 (j> can be calculated directly from the 
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definition of a timelike circular geodesic parameterized 
by the proper time s, leading finally to the following ex- 
pression for the worldline of a timelike circular equatorial 
geodesic in Schwarzschild coordinates 



x a (s) 



■ tt/2, 



M/r 3 



3M/r 



(6) 



Clearly, since the trajectory ^ is a geodesic, its acceler- 
ation a b = x c V c x h calculated with respect to the metric 
vanishes identically. However, we will proceed here in 
a different manner. In the next section, we will consider 
quantum effects as realized by observers with circular 
trajectories as ([6]) in the Minkowski spacetime. An ob- 
server following the worldline ([5]) in Minkowski spacetime 
experiences a centripetal acceleration 



M/r 



1 - 3M/r 



(7) 



On the other hand, in Minkowski spacetime the worldline 
of a particle in a uniform circular motion on the equato- 
rial plane with angular velocity Q is given by 

x a (§) = (t,r,ir/2,Clt), (8) 

from which one has immediately x a — 7(1,0,0, fi) and 



= v 



aba" 



rj 2 il 2 , where the constant 7 = dt/ds 



(1 — r 2 i7 2 ) -1 / 2 corresponds to the Lorentz factor. The 
angular velocity f2 is to be determined by imposing the 
centripetal acceleration ([7]) for the trajectory (JSJ) , yielding 



n = 



M/r 3 



y 1 - 2M/r 
and the following Lorentz factor 



7 : 



1 - 2M/r 
1 - 3M/r ' 



(9) 
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The treatment of the quantum effects realized by ob- 
servers in the circular geodesic (jSJ) of Schwarzchild space- 
time by means of an effective Minkowskian circular tra- 
jectory is, of course, only an approximation. It is shown 
in [12| , nevertheless, that the results obtained in a semi- 
classical approach assuming a Schwarzschild spacetime 
and a flat spacetime with external "Newtonian" attrac- 
tion forces such that ||7J) and ([9]) hold differ by no more 
than 30%, if we restrict ourselves to the circular orbits 
with r > 3M. Since Schwarzschild spacetime is asymp- 
totically flat, it is indeed natural that the emitted pow- 
ers calculated in Minkowski and Schwarzschild spacetime 
agree when one considers circular motions with large r. 
In fact, as it is shown in p[ and [lH, the power emitted by 
a particle in circular motion with radius r in Minkowski 
spacetime with angular velocity ([9]) is very close to that 
one emitted by a particle in a circular geodesic with 
the same radius r in a Schwarzschild spacetime, pro- 
vided that r > 6AI. This is in agreement with the well 



known fact that processes involving wavelengths with the 
same order of magnitude of the Schwarzschild radius need 
necessarily to be analyzed using fully curved spacetime 
calculations. Moreover, the acceleration defined in (JT]) 
has the additional desirable feature of being divergent at 
r = 3M, in accordance with previous works on geodesic 
emission which established that near the photonic 
orbit the emitted power diverges. The angular velocity 
([9]) mimics the main qualitative properties of the real 
Schwarzschild circular geodesies, justifying our assump- 
tion of circular trajectories in a flat spacetime with cen- 
tripetal acceleration ([7]). 



III. EMISSION RATES AND LIFETIMES 

The semiclassical current formalism employed in [J] to 
the proton decay case consists basically in considering the 
proton and the neutron as distinct energy eigenstates \p) 
and |n) of a two-level system such that H \p) = m p \p) 
and Hq\ti) = m n \n), where Hq is the proper Hamilto- 



nian of the system, 



and m n are, respectively, the 



proton and neutron masses. The weak channel ([I]) is 
implemented by considering a vector current associated 
to the two-level system coupled to a quantized fermionic 
field (corresponding to the electron e~ and to the anti- 
neutrino u e ) by means of the effective coupling constant 
G w , which is about the order of the Fermi coupling con- 
stant Gf ~ 1.166 x 10 _5 GeV -2 , whereas the strong 
channel @ involves a scalar current coupled to a quan- 
tized bosonic field (the pion ir~) by means of the effective 
coupling constant G s , of the order of the pion-nucleon- 
nucleon strong coupling g 2 NN /4Tr « 14[l4j. The currents 
are then specialized to the case of uniform circular tra- 
jectories with radius r and angular velocity f2 (and cen- 
tripetal acceleration a — rj 2 tt 2 ) in Minkowski spacetime. 

The proper decay rates corresponding to the weak ([1]) 
and to the strong @ channels are given, respectively, by 



G 2 a b 
' 8tt 4 



dXe is A{ j c)Zb T x 
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16 



l 2 Z a Z- 



and 



n — >p 



G 2 a 
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(11) 
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where n 2 = (m 2 +m 2 )/a 2 



and m v being, respectively, 



the electron and neutrino masses, and S = (m n — m p )/a. 
We assume here m„ = 0. The rates (fTTj) and (fT2|) are ob- 
tained, respectively, from equations (3.31) and (3.16) of 
, where all the relevant details can be found. In partic- 
ular, we have Z a = (-X+ie, 0, -(2ra/ 7 ) sin(OA7/2a), 0), 
where < e <C 1 is a regulator. For the relativistic case 
(7> 1), the terms involved in the integrations (|11[) and 
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HI are [I 



(13) 



Z a Z a » — L (A + iV3A + ^j (x + iV3A_ 



A - iV3B + ) (A - ^5+) , (14) 



where A± = l±yjl + 2e/ v / 3, B± = 1 ± y 1 - 2e/\/3, 
see Fig. [TJ 




FIG. 1: Path used in the complex integrations and (|12|l . 
Both terms (|17[) and (| 18j) come from the pole B+. They co- 
incide with the values calculated in [J by integrating around 
the pole A+. Notice that, for the processes considered there, 
the term 5 has a different sign. The second terms in (|15[1 and 
(|16|) come from the (degenerated) poles A- and B-. 

In contrast to the case considered in since m n > 
m p , one needs here to perform the complex integrations 
([IT]) and fljjj! along the path depicted in Fig. [TJ We get, 
after taking the limit e — > 0, 



- n — yp 



and 



-pw 

p — >n 



G 



2 „5 ; 



90tt 3 



(20 + 155 2 + 3<5 4 



2tt 



15^ 2 (l + <5 2 )) 
(15) 
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where Tp_^ n and r p _, n correspond, respectively, to the 
proper decay rates associated to the inverse processes 
considered in Hi, 



pW 

p — >n 



and 



G 2 a 5 



-2 v / 3<5 



49v / 3 + 102<5 + 30V3<5 2 



1728tt 3 

-u 2 (S9V3 + 905 + 36\/3(5 : 



12<5 3 
(17) 



2„ -2\/3<5 



/(8V37T) 



The approximations involved in the derivation of the ex- 
pressions pT|) and (|12p require, respectively, that the 



centripetal acceleration a obeys m e < a < m p and 
to,,. < a < TO p , where m, stands for the pion 7r~ mass. 
We notice that, for accelerations a ^> m p , the no-recoil 
hypothesis is violated^ and, hence, our approximation 
breaks down. The neutron proper lifetime associated 
with the decay rates (fTS"]) and ([16]) are given simply by 
= l/r™_^ p and = 1/r^ . Our results assume a 
particularly simple form if one considers the proton and 
neutron lifetime ratio, namely 



T%K = 1 



12<5e 



2V3i5 



and 



'pin 



£6 2^3* £(M 
5 Q(M)' 



(19) 



(20) 



where P(<$, ix) and Q(<5, /x) are polynomials easily obtained 
from (|15p and (fT7|) . It is clear that for large a, both ratios 
obey 



T p /T n ^l + 0(a 



(21) 



The strong channel is expected to be the dominant decay 
mode for neutrons with centripetal accelerations a such 
that < a < nip. Taking into account that m w ps 
139.57 MeV, m n ps 939.56 MeV and that m p ps 938.27 
MeV, we have that the proton and the neutron lifetimes 
differ by no more than 1% to 10% in the range of accel- 
erations where the strong channel dominates. The weak 
channel, on the other hand, dominates for smaller accel- 
erations m e < a < m^. (We remind that m e ps 0.51 
MeV.) From and it is clear that for 



a 3> a c — m n — m p ps 1.29 Mev 



(22) 



the asymptotic expression (|2T[) holds accurately. The 
meaning of the "critical" acceleration a c will be discussed 
in the last section. Here, we mention only that a c belongs 
to the range where the weak channel dominates. In fact, 
for m e < a < a c , a significative difference between the 
proton and the neutron lifetime is observed. 

A proper acceleration of the order of 1 MeV is ex- 
tremely high. For sake of comparison, protons in the 
CERN Large Hadron Collider have proper acceleration 
a ps 10~ 8 MeV[4j. In order to compare a c with cen- 
tripetal accelerations induced by realistic black- holes, we 
cast (|7J| in the form 



1.34 x 10 



-16 



Mo 
M 



(GM/rc 2 ) 
1 - 3GM/r 



(23) 



where a is now measured in MeV and the numerical con- 
stant corresponds to He 3 /GMq. Hence, realistic black 
holes (M ps Mq) will induce centripetal accelerations of 
the MeV order only for those circular orbits very close 
ng\ to the photonic orbit r = 3GM/c . Smaller black holes, 



nevertheless, can induce considerably higher centripetal 
accelerations for the (unstable) circular orbits located be- 
tween the photonic orbit and the last stable circular orbit 
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FIG. 2: Centripetal acceleration, in MeV, given by equa- 
tion (|23l) . for circular trajectories such that 3GM/c 2 < r < 
6GM/c 2 . The curves (a) to (d) correspond, respectively, to 
the ratios M/M Q = 10" 15 , 1(T 16 , 10~ 17 , and 10" 18 . We recall 
that, for the neutron decay, the approximations involved in 
the derivation of our results require a < m p ~ 938.27 MeV 
(no-recoil hypothesis) and a > m e ~ 0.51 MeV (weak chan- 
nel) or a > tUtt w 139.57 MeV (strong channel). The hypoth- 
esis of an extremely relativistic motion (7 ^> 1), on the other 
hand, requires orbits close to the photonic one, see equation 
(|24[). In the detail, the hatched area denotes the region of 
validity of our approximations, assuming 7 > 10. 



at r = 6GM/c 2 , see Figure [5] On the other hand, the 
hypothesis of an extremely relativistic motion (7 > 1) 
adopted in the approximations (| 13[) and (|14[) requires cir- 
cular trajectories close to the photonic orbit, since from 
(JTUJ) one has 



FIG. 3: Neutron proper lifetime r, in seconds, for circular 
trajectories such that 3GM/c 2 < r < 6GM/c 2 . The curves 
(a) to (d) correspond, respectively, to the ratios M/Mq = 
10" 16 , 10 -17 , 10~ 18 , and 10~ 19 . The (free) neutron inertial 
lifetime is approximately 886 s. The validity region of our 
approximations is depicted in Fig. 2. 
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Figure [3] depicts the lifetime for a neutron due to weak 
decay in circular orbits with 3GM/c 2 < r < 6GM/c 2 
around small black holes. As expected, the smaller is 
the black hole, the larger is the reduction in the particle 
lifetime. The semiclassical approach used here can be 
applied for other unstable particles as well. Particularly 
interesting is the muon weak decay pT — — > e^v e v^ which 
can also be described by a vector current based on a two- 
level system coupled to quantized fermions (the neutrinos 
v e and Uy) by means of a coupling constant of the same 
order than Gf- Since neutrinos have very small masses, 
the rate (fT"5)) for the muon weak decay is accurate for 
accelerations in the range < a < m e w 0.51 MeV. 
(We remind that m M w 105.7 MeV.) Figure 2] depicts 
the lifetime of a muon in geodesic circular orbits of small 
black holes such that ZGMjc 2 <r< 6GM/c 2 



FIG. 4: Muon proper lifetime r, in seconds, for circular tra- 
jectories such that 3GM/c 2 < r < 6GM/c 2 . The curves 
(a) to (d) correspond, respectively, to the ratios M/Mq = 



irr 



and 10" 



The muon inertial lifetime is 



about 2.2 x 10 s and the branching ratio corresponding to 
the inertial process fi~ — > e _ ^ e i^ M is greater than 98%. The 
validity region of our approximations is depicted in Fig. 2. 



IV. DISCUSSION 

As Figures [3] and [4] show, small black holes are neces- 
sary in order to induce sensitive alterations in the proper 
lifetime of unstable particles in circular orbits. However, 
noticeable effects do occur in the vicinity of the photonic 
orbit for realistic black hole. Despite that the analy- 
sis presented here is restricted to the (unstable) circu- 
lar geodesies close to r = 3M, it can give some hints 
about the behavior of particles in more realistic situa- 
tions. As equations (fT§|) and (p?0|) reveal, free protons 
and neutrons in geodesic circular motion close to the 
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photonic orbits have comparable proper lifetime. This 
situation is completely different from the inertial one, 
and its implication to particle physics in the vicinity of 
black holes has not been sufficiently studied yet. A sim- 
ilar conclusion holds for the muon. For a accelerations 
a~^> a c — 2\/3(«V — m e ) = 364.4 MeV, neglecting possi- 
ble effects of back- reaction]^, the muon and the electron, 
which in such a case can indeed decay by the inverse pro- 
cess e~ — > [i^v^v e , have comparable lifetimes. 

In order to grasp the meaning of the temperature ([3]), 
let us consider the case of the weak decay of protons 
and neutron in uniformly accelerated motion, where the 
Unruh temperature [To| Tjj — a/2n is known to play a 
central role[3|, Q . The proton and neutron lifetime ratio 
for this case can be obtained from the decay rates (3.13) 
and (3.17) of reference 0, 

T ;/ T w = (25) 

which also has the asymptotic form (|2"Tj) for large values 
of a. Notice that, in this case, we have exactly the same 
critical acceleration a c of (f2"2")l . 

In the case of the uniformly accelerated motion, one 
can describe the decay of protons and neutrons as seen 
by comoving Rindler observers. The key point here is 
that Rindler observers realize the inertial vacuum as a 
thermal state with temperature Tjj — a/2ir. Heuristi- 
cally, one can imagine the two-level system in thermal 
equilibrium with the Unruh radiation associated with the 
quantized fields in question. For our two-level system in 
thermal equilibrium at temperature T, the probability of 
occupation of the proton \p) and neutron \n) states are, 
respectively 

e -m p /T e -m n /T 
N P ~~ e -m p /T _|_ e -m„/T' Nn = e ~rn p /T + g -m„/T ' 

(26) 

The ratio N p /N n = e (™"- m p)/ T diverges for T -> 0, in- 
dicating that for low temperatures, the system is likely 
to be in its fundamental state. However, for tempera- 
tures T ^> (m n — m p ) the ratio tends to 1, indicating 
that the system can be in the states \p) or \n) with equal 



probability. In other words, the transitions \p) — » \n) 
and \n) — > \p) become equally probable for high temper- 
atures, shedding some light in the expression (|25| . For 
linear accelerations a such that the associate Unruh tem- 
perature Tjj — a/2n is much higher than the energy gap 
m n — m p , it is natural to expect that protons and neu- 
trons have the same lifetime, since both transitions of the 
two level systems are equally probable. The lifetime ra- 
tios (fT9|) and ([20]) suggest something similar to the case 
of uniform circular trajectories. They can be understood 
if one considers the two-level system in equilibrium with 
the quasi-thermal radiation with temperature © associ- 
ated with the quantized fields in question, confirming the 
view that observers in relativistic circular motion with 
centripetal acceleration a do realize the inertial vacuum 
as a quasi-thermal state with temperature ([3]) for the ex- 
tremely relativistic case[8, 9]. We stress that this state is 
quasi-thermal in the sense that it can be described by a 
temperature that varies slowly with the energy gap AE 
of the two-level system, monotonically from (n/2y3)Tjj 
to (n / y/3)Tjj , corresponding resp ectively to low and to 
high values of AE/a 0, [la, 13 • For the neutron de- 
cay, for instance, AE = m n — m p w 1.29 MeV, implying 
that for circular trajectories close to the photonic orbit 
r = 'iGM/c 2 , the quasi-thermal state is indeed charac- 
terized by a temperature close to the lower bound of ([3]). 
For the muon decay, on the other hand, the temperature 
is closer to the upper bound of (J3j) . A definitive answer 
to this problem, however, must necessarily face the sub- 
tle issue of quantum thermal distributions for rotating 
systems [13 ■ 



Acknowledgments 

The authors are grateful to G. Matsas for enlighten- 
ing discussions. This work was supported by FAPESP, 
CNPq, and UFABC. A.S. is grateful to Prof. V. 
Mukhanov for the warm hospitality at the Ludwig- 
Maximilians-Universitat, Munich, where part of this 
work was carried out. 



[1] R. Muller, Phys. Rev. D 56, 953 (1997); 
[2] G.E.A. Matsas and D.A.T. Vanzella, Phys. Rev. D 59, 
094004 (1999). 

[3] D.A.T. Vanzella and G.E.A. Matsas, Phys. Rev. Lett. 

87, 151301 (2001); Phys. Rev. D 63, 014010 (2000); H. 

Suzuki and K. Yamada, Phys. Rev. D 67, 065002 (2003). 
[4] D. Fregolente, G.E.A. Matsas, and D.A.T. Vanzella, 

Phys. Rev. D74, 045032 (2006). 
[5] V.L. Ginzburg and G.F. Zharkov, Sov. Phys. JETP 20, 

1525 (1965). 

[6] G.F. Zharkov, Sov. J. Nucl. Phys. 1, 120 (1965) 
[7] V.I. Ritus, J. Sov. Laser Res. 6, 497 (1985). 
[8] S.K. Kim, K.S. Soh, and J.H. Yee, Phys. Rev. D 35, 557 
(1987). For earlier works on the subject, see: JR. Letaw 



and J.D. Pfautsch, Phys. Rev. D 22, 1345 (1980); JR. 

Letaw, Phys. Rev. D 23, 1709 (1981); JR. Letaw and 

J.D. Pfautsch, Phys. Rev. D 24, 1491 (1981); J. Math. 

Phys. 23, 425 (1982). 
[9] L.C.B. Crispino, A. Higuchi, and G.E.A. Matsas, The 

Unruh effect and its applications, to appear in Rev. Mod. 

Phys. |arXiv:0710.5373] . 
[10] W. G. Unruh, Phys. Rev. D 14, 870 (1976). 
[11] R.M. Wald, General Relativity, University Of Chicago 

Press (1984). 

[12] L.C.B. Crispino, A. Higuchi, and G.E.A. Matsas, Class. 
Quantum Grav. 17, 19 (2000); J. Castineiras, L.C.B. 
Crispino, G.E.A. Matsas, and R. Murta, Phys. Rev. D 
71, 104013 (2005). 



6 



[13] C. W. Misner, R. A. Breuer, D. R. Brill, P. L. 
Chrzanowski, H. G. Hughes III and C. M. Pereira, 
Phys. Rev. Lett. 28, 998 (1972). R. A. Breuer, P. L. 
Chrzanowski, H. G. Hughes III and C. W. Misner, Phys. 
Rev. D 8, 4309 (1973). R. A. Breuer, Gravitational Per- 
turbation Theory and Synchrotron Radiation - Lecture 
Notes in Physics (Springer-Verlag, Heidelberg, 1975). 



[14] R. Machleidt and I. Slaus, J. Phys. G 27, R69 (2001). 
[15] W. G. Unruh, Phys. Rep. 307, 163 (1998). 
[16] N. Obadia and M. Milgrom, Phys. Rev. D 75, 065006 
(2007). 

[17] G. Duffy and A. C. Ottewill, Phys. Rev. D 67, 044002 
(2003). 



